Soliton dynamics in coupled system of scalar fields with broken U(l) symmetry 
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We construct a C/(l)-breaking nonlinear system which bears topological solitons with in- 
■ teresting properties. The symmetry is broken both spontaneously and exphcitly. The model 

■ 

presented has TWO degenerate pairs of kinks and antikinks which can be transformed into 



each other by perturbations induced via external triggers and/or collisions. We calculate 
both the single soliton solutions and collision processes. These include simple scattering, 
pair annihilation, pair capture, and excitation-de-excitation processes. 
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I. INTRODUCTION 

Q^ ■ There are some noticeable similarities between relativistic solitons and classical particles. For 

l/^ ' instance, interacting with each other via short range forces and collide without losing their identities 

^ ■ firi 

^sj ' llrl5l|) ^-Iso they are localized and do not disperse while propagating in medium. Because of their 

irl ' field nature, they do tunnel a barrier in certain cases, although this tunnelling is different from the 

well-known quantum version [3, |6|. Topological solitons are stable, due to the boundary conditions 
at spatial infinity. Their existence, therefore, is essentially dependent on the presence of degenerate 
vacua 

Topology provides an elegant way of classifying solitons in various sectors according to the 
mappings between the degenerate vacua of the field and the points at spatial infinity. For the Sine- 
Gordon system in 1 -|- 1 dimensions, these mappings are between (p = 2mr, n G Z and x = ±00, 
which correspond to kinks and antikinks of the SG system. More complicated mappings occur 
in solitons in higher dimensions Q|. Coupled systems of scalar fields with soliton solutions have 



found interesting applications in double-strand, long molecules like the DNA molecules 
bi-dimensional QCD 1^ . 
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Systems of two coupled real scalar fields with particular self-interaction potentials were inves- 
tigated by Bazeia et al 3| and Riazi et al 0, [l^ . Inspired by the coupled system introduced in 
0,13] J we propose a new coupled system of two real scalar fields which shows interesting types of 
solitons with well-defined topological charges and rest energies (masses). The present model may 
be applied to the dynamics of domain walls in the hybrid inflationary cosmology [igI]. 

In this paper, we focus on a system with the double scalar field potential 

1 



(1) 



in which (j) and ip are real scalar fields, and A is a constant and the explicit symmetry breaking 
parameter as well. This potential which is similar, but not quite the same as that of the hybrid 
inflationary model[l^, is plotted in Fig. [T] for A = 1. Note that the potential along the cj) axis 
has always two degenerate vacua at (j) = if) while the potential along the ip axis has minima at 
ip = ±1 only if A < 4. For A > 4, the potential has only an absolute minimum at ip = (see Figure 

ED. 




FIG. 1: The self-interaction potential (1) is shown as a contour map on the {(pjip) plane. The two bluish 
sites indicate the degenerate minima of the potential. 



II. LAGRANGIAN, FIELD EQUATIONS AND SYMMETRIES 



Within a relativistic formulation, the Lagrangian density of our proposed system is given by: 

1 



(2) 



FIG. 2: Plots showing V{(f), 0) and V{0, V') for A = 0, 1, 10 (from bottom to top). 

in which A is the U (1) exphcit symmetry breaking parameter. One can write this Lagrangian in 
terms of the complex scalar field <I>, where 

^ = (l) + iip, (3) 
In terms of this complex scalar field the Lagrangian density of system reads 

C = la'^^t^^^ _ ($t$ _ 1)2 ^ Ixilm^f. (4) 

From either of these two forms of the Lagrangian density, the following field equations for (j) and 
ip are obtained: 

□ (/.= -40((/.2+V'2-l); (5) 

and 

Dip = + - 1) + Xip. (6) 

It is obvious that if A = 0, the Lagrangian density (l2|) or is both Lorentz invariant and also 
invariant under a global U{1) symmetry 

$ ^ = (7) 



or 



4>' = (j) cos 9 — Ip sin 0, 

Ip' = (psm. e + cos 6. (8) 
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The corresponding energy-momentum tensor[17l. Il8l| of the system is: 

Tf,u = df,4>d,y(t) + d^ijdui} - g^,vC; (9) 
which satisfies the conservation law 

a^T^'^ = 0. (10) 

In Equation ([9]), g^u = diag{l, —1) is the metric of the (1 + l)-dimensional spacetime. The Hamil- 
tonian (energy) density is obtained from Eq. (jlip according to 

Like other systems bearing topological solitons, the present system also has the following topological 
current: 

JM = ^e^-d,cP (12) 

which is locally coserved: 

df.J" = 0. (13) 
The corresponding topological charge is given by: 

r-+oo 1 

Q= / jOrfx = -[0(+oo) -0(-oo)]. (14) 

J ~oo ^ 

Note that since the vacua of the system reside at {4>,ip) = (±1,0), only the 0- field is responsible 
for the topological charge. 

According to the Goldstone theorem, if a continuous global symmetry is broken spontaneously, 



IS 



20l |. However, 



there appears a massless (Goldstone) particle for each broken group parameter 
in the case of the Lagrangian density Q, in addition to the spontaneous breaking of the U(l) 
symmetry for A = 0, the symmetry is broken explicitly by the A-term. If we expand the potential 
around either of the vacua (0 = ±l,'i/' = 0), there appears the following mass terms: 

V{x, V") - ^m^X^ + ]^m^i?, (15) 

where 



2 l(</'=±i,v=o) - o> 



(16) 
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and 



^i' = "^l(</'=±i-V'=o) = (17) 



where x = 4'~^- It is seen that due to the exphcit symmetry breaking term, the massless Goldstone 
boson which is normally massless, has acquired a mass (A). 

Let us consider the C/(l)-symmetric case A = 0. The potential reduces to the well-known 
complex (f"^ model and we have the global U{1) symmetry 

$ — > $ = e*^$, 

C = C. (18) 

This symmetry leads to the following conserved current and charge, as deduced from the celebrated 
Noether's theorem[l8] 

= i{^*d''^-^d^^*). (19) 



Qn = I J%dx. (20) 



Writing the complex scalar field ^ in the form <^> = i2e*^, the current p9p and the charge ()20p take 
the following simple forms: 

= 2R^d^i, (21) 

and 

/ + 00 
R^d^idx. (22) 
-oo 

It is obvious that the Noether charge vanishes for all static solutions, including the static kinks 
and antikinks to be introduced shortly. For a time-varying field like $ = -R(x)e*'^*, however, we 
have the non- vanishing Noether charge Qn = 2ijj J R?'{x)dx. 

By using the dynamical equations, it can be easily shown that the U{1) current is partially 
conserved and we have 

d^r^ = 2A<^V, (23) 

which is proportional to A, like the situation arising in PCAC (=partially conserved axial current), 
although in PCAC the explicit symmetry breaking term in the Lagrangian is usually assumed to 
be linear in T^jisj]. 
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III. KINK AND ANTIKINK SOLUTIONS 

In order to find the static kink and antikink solutions, we have employed the following numerical 
procedure within the MATLAB environment. The program starts with an approximate solution 
which is the exact solution of the ?7(l)-symmetric system (A = 0). The solution is then perturbed 
via small changes in the field values and the total energy of the system, as given by 

E = j Hdx (24) 

is calculated at each step. The small changes in the field values is accepted if the total energy 
is reduced in each step, otherwise it is rejected. This procedure is iterated repeatedly, until the 
program reaches a minimum energy configuration. Kink and anti-kink solutions obtained in this 
way are shown in Figures [3l The corresponding orbits on the ((/>, ^)-plane are shown in Figure [H 
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FIG. 3: Minimum energy, static solutions. The solid curves represent <^ and the dash-dotted curves are for 



IV. BOGOMOLNY CONDITION 



It is well known that in many nonlinear equations bearing topological solitons, static solutions 
satisfy the so-called Bogomolny condition. This condition puts a lower bound on the total energy 
of the system which is proportional to the topological charge 
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22 . 



In this section, the Bogomolny bound is calculated and its proportionality to the topological 
charge is checked. 
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FIG. 4: The minimum energy, static solutions mapped on the {4>,il') plane. 



Multiplying equation ([5]) by cp' and adding it to equation ([6]) multiplied by ^p', we obtain 



(25) 



Here, prime means derivative with respect to x. We thus obtain the following first integral of the 
static field equations: 



(26) 



where C is a constant of integration. For localized field equations, this constant should be zero, 
since the fields rest on their vacuum values at x — t- ±00. 

The above result can also be obtained from the following Bogomolny-type argument. One starts 
with 

2 



1 1 \ 

^(/>' + ^V'±V^ >o. 



Squaring and integrating over x, we obtain 



+00 / 1 1 



Introducing a superpotential "0) such that 



(<^i,V')^ 



/2y ((i(/> + di\)) 



we obtain 



dW dW , — 
dW = —d^ + —di^ = V2V{d(t> + #), 
0(p oip 



E > \W+-W-\, 



(27) 



(28) 



(29) 



(30) 



where W± = W {(j){zizoo) , 11^(^00)) . The Bogomolny bound is saturated for the equality sign in (|27p . 
Since the vacua of the system under consideration reside on = ±1 and -0 = 0, it is the (j) field 
which is responsible for the topological charge. 
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V. INTERACTION BETWEEN KINKS AND ANTIKINKS 



As in other kink-bearing systems, an important question is the form of the inter-kink potential 
and the behavior of the sohtons in cohisions with each other. Kinks and antikinks of different 
nonlinear systems behave differently in collisions. In most cases, the following situations arise: 1) 
A pair of kinks or antikinks which have the same topological charge repel each other. They retain 
their original shape after the collision. In the sine-Gordon system which is integrable, the pair 
retain their original speeds after the collision. In non-integrable systems like the (j)'^ system, part 
of the energy is converted into small amplitude waves which are radiated away and the final speed 
of the solitons is less than their initial speed. 2) In the sine-Gordon system, the collision between a 
kink and an antikink does not lead to their destruction and the pair retain their initial speeds after 
the collision. The force between the pair is velocity-dependent j^]. It is attractive at relatively 
large distances and repulsive at short distances. In non-integrable systems like (p'^ or double-sine- 
Gordon system, the collision process between a kink and an antikink is more complicated and 
interesting phenomena happen. For example, the pair annihilate each other when their relative 
velocity is smaller than a first threshold vi . For velocities larger than vi and smaller than a second 
threshold V2, there appear scattering windows in which the pair leave the interaction region with 
a smaller speed than their initial speed. Some small amplitude waves are radiated away in this 
process. Velocities larger than V2 lead to the scattering of the pair and emission of radiation. 3) In 
the sine-Gordon system, there is a bound state solution in which a kink and an antikink oscillate 
around the center of mass of the system indefinitely. Breather solutions in non-integrable systems 
like 4>'^ are unstable and lead to the annihilation of the pair after transient oscillations. 

In this section, we are interested in studying the above possibilities in the system under investi- 
gation. Since analytical calculations are not possible here, we employ the finite-difference method 
described in {2]. 

Figure [5] shows simple scattering of a kink and antikink. The velocity of each soliton (in units of 
the speed of light) is 0.6 for this process. Figures [6][8] show examples of some interesting interactions 
for the system considered in this paper. In Figure [6l the pair annihilate each other into a pair of 
neutral wave packets which leave the interaction area with larger velocities. Figure [7] shows the 
formation of a bound pair (or soliton molecule), emitting the residual energy in the form of lower 
amplitude waves. In Figure [HI the collision leads to the excitation of each soliton and the excited 
soliton relaxes into its lower energy state, emitting a neutral wave packet within a short time. In 
these plots, positive topological charge density is shown in red and negative charge in blue in order 
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to better illustrate the fate of charged objects. 
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FIG. 5: Kink-antikink simple scattering at u = 0.6. Topological charge density is plotted on the {x, t) plane. 
Reddish color indicates positive and bluish indicates negative topological charge. 
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FIG. 6: Annihilation of a kink-antikink pair into a pair of neutral wave packets at v = 0.36. Reddish color 
indicates positive and bluish indicates negative topological charge. 

Let us consider the soliton molecule in more detail. If there is a stable kink-antikink bound 
state, then one should be able to obtain it via an energy-minimization procedure. To this end, 
we have written an energy-minimization program in the MATLAB environment, which produces 
a minimum-energy solution, starting with a trial pair of functions which satisfy the boundary 
conditions and the general functional form of the soliton molecule. The initial guess functions read 

m = - 1, (31) 



and 



^(^) = -rr-2- (32) 



10 



1000 2000 



4000 5000 6000 

Time Steps 



7000 8000 



FIG. 7: Formation of a soliton molecule via kink-antikink capture at v — 0.5. Note that the surplus energy 
is radiated away. Reddish color indicates positive and bluish indicates negative topological charge. 
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FIG. 8: Kink-antikink interaction, leading to the excitation and subsequent de-excitation and recoil of the 
kink and antikink {v — 0.7). Reddish color indicates positive and bluish indicates negative topological 
charge. 

It is obvious that this trial equation has zero total topological charge, comprising equal negative 
and positive charge densities of the kink and antikink constituents. The initial guess, together with 
the minimum energy solution are shown in Figure [9j This minimum-energy bound state of the 
kink and antikink closely resemble the pair formed in the numerical experiment shown in Figure El 



VI. CONCLUSION 



Kink-bearirig systems show very interesting phenomena both in numerical calculations and in 
the real world [iHgI]. A complex scalar field with U{1) symmetry is well known and worked out 
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FIG. 9: Initial guess functions ([31]) and ([5^ shown as dashed curves, together with the minimum energy 
kink-antikink bound state solution (solid lines). 



thoroughly in field theory 18t|. When the global U{1) symmetry is made local, it leads to the 
appearance of electric charge and electromagnetic interactions. Motivated by these interesting 
properties, we constructed a double-real-field Lagrangian with a C/(l)-breaking term. We obtained 
static kink and antikink solutions and showed that there are two degenerate pairs of kinks and 
antikinks in the system. Several numerical experiments were performed to explore what happens in 
the kink-antikink collision at various relative velocities. It was observed that different interesting 
phenomena may happen. Examples include simple scattering, pair annihilation into neutral wave 
packets, formation of soliton molecule (bound kink-antikink pair), and excitation-decay process. 
The soliton molecule formed in some kink-antikink collisions approximately conforms with the 
solution obtained via minimizing the energy of a pair of guess functions adapted to the required 
topological charge and boundary conditions. In order to distinguish charged solitons from neutral 
wave packets and follow the evolution of each charged soliton, we preferred to plot charge density 
rather than the fields or energy densities which is usually practiced in the literature. 

One final observation to be pointed out is that in many examples, the resulting dynamics is not 
symmetrical about the pair center of mass. In other words, there is a left-right asymmetry which 
constitute yet another difference with other well-known non-integrable systems. 

The results presented in this paper could probably be applied to the dynamics of domain walls 
in hybrid inflationary models with modified potentials 1^ . 
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